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The one-loop nuclear structure corrections of order (Za)^ to the Lamb shift and hyperfine 
splitting of the deuterium are calculated. The contribution of the deuteron structure effects to 
the isotope shift (ep)-(ed), (pp) — (pd) in the interval IS -i- 2S is obtained on the basis of modern 
experimental data on the deuteron electromagnetic form factors. The comparison with the similar 
contributions to the Lamb shift for the electronic and muonic hydrogen shows, that the relative 
contribution due to the nucleus structure increases when passing from the hydrogen to the deuterium. 


I. INTRODUCTION 

Experimental and theoretical investigation of the Lamb shift and hyperfine structure of muonic hydrogen (pp) and 
muonic deuterium (/rd) can lead to essential progress in determining important fundamental parameters of the proton 
and deuteron. A better understanding of the nuclear structure and polarizability effects in these hydrogenic atoms 
can be gained due to such studies [1-5]. Namely the proton structure and polarizability corrections lead to main 
theoretical uncertainties on the Lamb shift and in the hyperfine splittings (HFS) of the hydrogen atom. One of the 
recent studies of these effects in the HFS in electronic and muonic hydrogen carried out in [6] shows that 70 % of the 
difference (AE^™ — AE^pg) = 0.046 MHz (QED contribution AE^p^ doesn’t take into account the proton recoil, 
structure and polarizability corrections) in the case of electronic hydrogen can be explained on the basis of effective 
field theory describing the interaction of baryons with photons and leptons. The measurement of the 2P -L 2S Lamb 
shift in muonic hydrogen has been carried out during recent years at PSI (Paul Scherrer Institute) [7]. The goal of 
the experiment is to measure the Lamb shift in muonic hydrogen with 30 ppm precision and to obtain the root mean 
square (rms) proton charge radius with 10“^ relative accuracy that is an order of the magnitude better than in the 
analysis of the elastic e-p scattering and the Lamb shift in electronic hydrogen. Another important task is connected 
with the study of the hydrogen - deuterium isotope shift for the interval IS -L 2S [8-10]. The experimental value of 
the H-D isotope shift 


AEH-D(lS-b2S) = 670 994 334.64(15) kHz, ,5 = 2.2 • 10-^° (1) 

was obtained with such high precision that the nucleus structure and polarizability effects should be taken into account 
on close theoretical examination of this quantity. Using the result (1) the numerical value for the difference of the 
proton and deuteron charge radii can be evaluated as follows [11]: 

r2-r^ = 3.8213(11.7) W. (2) 

So, when the Lamb shift measurement in the system /rp with the precision 30 ppm will be accomplished the deuteron 
radius rj can be derived by means of (2) with the accuracy 10“^. Another approach to determine more exact value 
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of the deuteron charge radius is related to the measurement of the muonic hydrogen - muonic deuterium isotope 
shift for the interval IS 2S. In this case the precision of the appropriate experiment should be comparable with 
(1). The possibility of such measurement of the /xp - /xd isotope shift for the interval IS2S is considered to be 
quite real taking into account the relative error of order 10“® at PSI for the Lamb shift measurement in muonic 
hydrogen. Evidently the calculation of all possible corrections to the Lamb shift in muonic hydrogen and deuterium 
with similar precision must be performed. The total value of the Lamb shift is determined by the sum of quantum 
electrodynamical (QED) contributions (one-loop, two-loop, three-loop corrections, recoil corrections, radiative recoil 
corrections) and nuclear corrections. Whereas QED contributions were obtained presently with the relative accuracy 
10“^ due to numerous calculations [11] the nuclear structure corrections are known less precisely. The accuracy of the 
corresponding calculation of the nucleus structure and polarizability effects depends on the experimentally measured 
nuclear densities of the charge, magnetic moment, ... and on the used theoretical models. Both experiments for the 
isotope shift (1) and the hydrogen HFS can play the selecting role among numerous nuclear models. 

The ground state hyperfine splitting in deuterium represents another important quantity where the nuclear structure 
corrections can be tested experimentally. The experimental value of the deuterium HFS was obtained with high 
accuracy many years ago [ 12 , 11 ]: 

AEHPg(D) = 327 384.352 521 9(17)(3) kHz, (5 = 5.2 x IQ-^^. (3) 

The difference between the experimental (3) and theoretical values for the deuterium HFS accounting for the QED 
corrections is equal to AE^pg(D) - AE^pg(D) = 45 kHz. The essential contribution to the theoretical quantity 
AE^ps(D) is given also by the nuclear structure corrections. The deuteron is the spin 1 particle and its electromagnetic 
structure is described by three form factors. The aim of our study consists in the exact consideration of the deuteron 
structure corrections of order (Za)® in the Lamb shift and HFS of the electronic and muonic deuterium. In distinction 
to previous investigations of these problems [13-16] we use the explicitly covariant approach to the description of the 
electromagnetic deuteron - lepton interaction to construct the quasipotential of the two-photon interaction. Moreover, 
we take into account modern experimental data on the deuteron electromagnetic form factors: charge monopole, charge 
quadrupole and magnetic dipole. We leave aside the polarizability effects mentioned above. 


II. CORRECTIONS OF ORDER {Zaf TO THE DEUTERIUM LAMB SHIFT 

The main nuclear structure dependent contribution of order (Za)^ to the hydrogen Lamb shift is determined by the 
one-photon interaction. In the one-photon exchange approximation the amplitude of the scattering process ed ^ ed 
is just the contraction of the electron and deuteron electromagnetic currents, multiplied by the photon propagator. 
The parameterization of the deuteron electromagnetic current takes the form [17,18]: 


Jd(P2,q2) = £*M2)< 


2m2 


(P2 + q2)/x kpk,, 


- 5:;;g^F3(k^) p,(P2). 


2m2 


2 m 2 


k^ 

^^" 2 m 2 


(4) 


where p 2 , q 2 are four momenta of the deuteron in the initial and final states, k = q 2 — P 2 , m 2 is the deuteron mass. 
The spin 1 polarization vectors satisfy the following conditions: 


A)£'"(k, A') =-(5 aa', k^e''(k. A) = 0 , ^ e* (k, A)ej.(k, A) =-g^^,-b 


k„k; 




(5) 


The generator of the infinitesimal Lorentz transformations 


( 6 ) 


The deuteron electromagnetic form factors Fi(k^) depend on the square of the photon four momentum. They are 
related to the deuteron charge Fc, magnetic Fm and quadrupole Fq form factors by the expressions: 

Fc = Fi -b -?7 [Fi -b (1 -b ? 7 )F 2 — F 3 ], 


Fm = F3, rj = — 


k" 

4m2 ’ 


(7) 
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Fq = Fi + (1 + ri)F2 — F 3 . 
The lepton electromagnetic current has the form: 


Jr(pi>qi) = u(qi) 




2 mi 


2 mi 


i(Pi), 


( 8 ) 


where pi,qi are the electron (muon) four momenta in the initial and final states, = ( 7 ^ 7 *^ — 7 '^ 7 ^)/ 2 , ki is the 
lepton anomalous magnetic moment, mi is the lepton mass. To obtain the one-photon interaction contribution to the 
Lamb shift, it is necessary to average the currents (4) and ( 8 ) over the electron and deuteron spins. As a result the 
contribution to the Lamb shift of order (Za)^ is expressed through the deuteron charge radius rd as follows: 


[4 + Fm( 0) - Fo(0)] , r» = 


(9) 


where /i = mim 2 /(mi -|- m 2 ) is the reduced mass. The numerical value of this expression for the interval IS q- 2S in 
muonic deuterium for the rd = 2.094 fm is equal to -186.74 meV. The magnetic quadrupole term of Eq. (9) which is 
proportional to the difference [Fm(0) — Fc(0)] = 2/id — 1, coincides with the result of Ref. [19] for the spin 1 particle. 
Its contribution to the energy spectrum amounts to 0.2 %. Consider the two-photon exchange amplitudes shown in 
fig.l. They give the corrections of order (Za)^ to the deuterium Lamb shift. The spin 1 particle (deuteron) exists in 
the intermediate state of the processes of the virtual Compton scattering. The amplitudes of the virtual Compton 
scattering on the lepton and deuteron can be presented as follows: 
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FIG. 1. Nuclear structure corrections of order (Zo)®. 


To extract the contribution of the two-photon amplitudes to the Lamb shift we can average the tensors (10) and 
(11) over the lepton and deuteron spins. Multiplying obtained relations and contracting over the Lorentz indices we 
can write necessary quasipotential in the form (we used the system FORM [20] for the trace calculations of the Dirac 
7 matrices and the contraction over the Lorentz indices): 
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[3m2 - + kg] + 4FiF2(k^ - kg)^[k^ - m 2 (k^ + 2kg)] + 4FiF3m2k^(k^ - kg)(k^ - 2kg)- 

[k^ - m2(k2 + Skg)] - 2FiF3k4(k2 - klf + Flk'^kl[-3k^ + 4m2(2k2 + kg)]!. 
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After rotating the kg contour of the Feynman loop integration we integrate Eq. (12) over the four-dimensional 
Euclidean space using the relation: 


d^k = 47r / k'^dk / sim 0 • d(/), ko = kcos0. 

Jo Jo 


(13) 


Then the integration over angle (j) in Eq.(12) can be done analytically. Calculating the matrix element of obtained 
potential between Coulomb wave functions we find the following one-dimensional integral representation for the 
contribution to the energy spectrum: 
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where hi = i/k^ -|- 4m?. The analysis of the coefficients with the deuteron form factors in Eq. (14) at small values of 
the variable k shows that there are three terms proportional to F^, F1F2, F1F3 which contain infrared divergences. 
For the regularization of these terms it is necessary to supplement the quasipotential of Eq. (14) by the iteration term 
which gives the following contribution to the energy levels: 

AE?(|, =< [Vlg, X X Vig,]^" >= (15) 
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where = (b^ — p^)/2/xr is the inverse free two-particle propagator [21], Vi-y is the quasipotential of the one- 

photon interaction. Furthermore to eliminate the divergence in the part ~ Fi(k^)F 3 (k^) we subtract from it the 
corresponding contribution accounting for the point-like deuteron ~ Fi( 0 )F 3 ( 0 ) (needless to say that the one-loop 
corrections of order (Za)® for the point deuteron in the Lamb shift must be studied independently [22], the relating 
studies are in progress). To perform the numerical calculations on the basis of Eq. (14) we employ recent experimental 
data on the deuteron electromagnetic form factors [23] extracted from the tensor polarization data in the elastic 
electron - deuteron scattering at given values of the four- momentum transfer. Some useful parameterizations exist 
for the deuteron form factors in the range 0 < k < 1.4 GeV. We used the parameterization from Ref. [24] of the 
following form: 



Fq I = Ff 


D 


M{ri) ■ 



where Fd is the nucleon dipole form factor, M{ri) is the matrix of the coefficients, ry = k^/dm^ and 


f„. = k“^ 


Q;^ • + 
[—1 '^mi ^ ^ 


(16) 


(17) 


The parameters ami, ami describing each electromagnetic form factor can be found in Ref. [24] (see also http://www- 
dapnia.cea.fr/Sphn/T20/Parametrisations). It is necessary to point out that the parameterization (16) leads to the 
following value of the deuteron charge radius rd = 2.094 fm [23] which is 2.6 % smaller than the value rd = 2.148 fm, 
obtained on the basis of Eq. (2) by means of the rp from Ref. [11]. The results of the numerical integration in Eq. 
(14) are presented in Table 1. The momenta of the integration which give the main contribution in Eq. (14) have 
a characteristic nuclear scale. The obtained correction of order (Za)^ 0.53 kHz for the 2S IS interval in electronic 
hydrogen is in agreement with the result 0.49 kHz of Ref. [16]. In the case of the muonic deuterium analogous 
contribution is equal to 2.57 meV. Then the value of the nuclear structure corrections to the isotope shift in muonic 
(electronic) hydrogen accounting the results of Ref. [25,26] is as follows: 


AE^^,(1S-2S) = 


/xp — fid 1.41 meV 
ep — ed 0.497 kHz 


(18) 


It must be taken into account when the comparison with the experimental data will become available. 


III. CORRECTIONS OF ORDER {Zaf TO THE DEUTERIUM HYPERFINE SPLITTING 

To construct the HFS part of the one-loop quasipotential one needs to keep not only the terms of the operator ( 6 ) 
with Sjj = 2 ieiji:S 2 which are proportional to the deuteron spin S 2 . The terms containing the operator Spn, which 
are expressed through the generator of the Lorentz boosts S for the spin 1 particle, also must be taken into account. 
There are two possibilities to determine the HFS part of the interaction operator: 

1. To keep consistently all the terms, which can contain the spin - spin interaction operator (S 1 S 2 ) when calculating 
the contraction of the amplitudes (10) and (11) over the Lorentz indices. 

2. To use the special projection operators on the states of the lepton and deuteron with the total spins 3/2 and 1/2. 
In this study we carry out covariant construction of the HFS quasipotential introducing the projection operators 
#^ 3/2 and 7 r^,i /2 for the particles in the initial and final states of the following form: 

TABLE I. Nuclear structure corrections of order (Za)"^, (Za)® to the Lamb shift in electronic and muonic deuterium. 




IS (Za)® 


2S (Za)® 



ed (kHz) 

6.875 • lO-" 

-0.603 


-0.075 


0.527 

pd (meV) 

213.42 

-2.94 

1 26.68 

-0.37 

1 -186.74 

2.57 
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where the spin-vector ^'/^(P) and the spinor d>(P) describe the lepton-deuteron states with the total spin 3/2 and 1/2 
respectively, M = mi -|- m 2 is the total mass, P = pi -I- P 2 - 

Multiplying amplitudes (10) and (11) and taking into consideration the relations (19)- (21) we can represent the 
necessary HFS quasipotential by the following expression: 
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The number of the form factor terms was reduced by half in comparison with Eq. (14) because the spin-dependent 
terms of the second particle are proportional to the form factor F 3 . The expression (22) is less singular than the 
operator (14). The sole infrared divergence of Eq. (22) when k ^ 0 is connected with the term ~ FiF 3 k^. It can be 
eliminated entirely by subtracting the iteration part of the quasipotential which can be obtained using Eqs. (4) and 
( 8 ): 
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Subtracting (23) from (22) we can make the analytical integration over the angle variables in the Euclidean momentum 
space. The averaging of the obtained expression over the Coulomb wave functions will lead to the appearance the 
factor | 7 /( 0 )p. As a result the contribution of the two-photon amplitudes to the deuterium HFS can be written as the 
one-dimensional integral: 
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where the Fermi energy of the ground state hyperfine splitting in deuterium 


Eq = 2fida 


4 M 


memp 


= 326 967.678(4) kHz, 


(25) 


fj,d= 0.857 438 228 4(94) is the deuteron magnetic moment in the nuclear magnetons, /j is the reduced mass of the 
deuterium atom, mp is the proton mass. The numerical integration of Eq. (24) was also done by means of the deuteron 
form factor parameterization (16). As a result the contribution of Eq. (24) to the HFS of the electronic deuterium 


EHFS(ed) = -34.72 kHz. 


(26) 
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The calculation of the deuteron structure corrections to the HFS of the ground state was performed previously in 
the analytical form in the zero radius approximation for the deuteron form factors. In this approximation the wave 
function of the deuteron D-state was omitted and the S-state wave function was written in the asymptotic form Be“'®“'. 
The contribution obtained in Ref. [15] is the following: 


AE^^®(ed) = • 
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a^(l + 21n2) 


3a me j K 

Stt mp me 



(27) 


3ck ixip , /cl/o i-\ \ 

- -In-(^2 _ 2^ 3 + ^2^ 

47r mr, vri ,,. \ p p / 
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-21.31 kHz. 


where k = 45.7 MeV is inverse deuteron size. The essential growth of the nuclear structure correction of order (Za)^ 
to the deuterium HES in our case as compared with Ref. [15] can be explained by removing the restriction of the zero 
radius approximation. The value E^.^® obtained here leads to significant growth of the difference between the theory 
and experiment which amounts to 45 kHz without consideration of the deuteron nuclear structure and polarizability 
contributions. In the case of muonic deuterium the corresponding contribution 

Ef.f®(/rd) = -0.925 meV. (28) 

The theoretical uncertainty of the results (18), (26), (28) and presented in Table 1 is determined by the errors of the 
experimental data for the deuteron electromagnetic form factors which amounts to 5 % in the most important interval 
0 < k < 0.5 Gev. So, the theoretical error of the results obtained in this work may amount to 10 %. It may be useful 
to compare the relative values of the corrections connected with the nuclear structure in the energy spectra of light 
and heavy hydrogen. The main one-loop contribution to the HFS of electronic (muonic) hydrogen is determined by 
the following expression (the Zemach correction) [27,28]: 


AE = Ep • Sz 



dp 

(p2 + b2)2 


Ge(-p^)Gm(-P^) 
1 -I- «: 


Ef(— 2 /xa)Rp, b = a/i. 


(29) 


where Rp is the Zemach radius. In the coordinate representation the Zemach correction is determined by the magnetic 
moment density pM(r) and by the charge density pE(r). The value Rp can be obtained in the numerical form by using 
the parameterization for the proton electromagnetic form factors obtained at Mainz 20 years ago from the analysis 
of elastic electron - proton scattering [29]. The relative contributions of the Zemach correction to the hydrogen HFS 
are as follows: 


electronic hydrogen : Rp = 1.067 fm, Sz = —40.3 ppm, (30) 

muonic hydrogen : Rp = 1.064 fm, Sz = —74.7 • 10“^. (31) 

The relative deuteron structure contributions obtained in this study have the following numerical values: 

electronic deuterium = —106.19 ppm (32) 

muonic deuterium = —188.37 • 10“^. (33) 


Three times growth of the contributions (32), (33) as compared with the expressions (30), (31) can be caused by 
increasing the distribution region of the deuteron electric charge and magnetic moment. The numerical results 
obtained in this work for the muonic deuterium must be taken into consideration when both extracting the deuteron 
charge radius in future experiments on the isotope shift (/ip) — (/id), and comparing theoretical predictions with 
measurements of the HFS in deuterium. 
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